Search for: All records

Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e. proto-galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: (i) can have a stellar mass range between 1 × 108 < M⋆ < 2 × 1010 [M⊙], a virial mass range between 3 × 1010 < M⋆ < 6 × 1011 [M⊙], and be as young as 8 ≲ Age ≲ 12.8 [Gyr] (1 ≲ z ≲ 6); (ii) are pre-dominantly centrally concentrated, with $\sim 50~{{\ \rm per\ cent}}$ of the stars contained within 5–10 kpc; (iii) on average show weak but systematic net rotation in the plane of the host’s disc at z = 0 (i.e. 0.25 ≲ 〈κ/κdisc〉 ≲ 0.8); (iv) present [α/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host’s old disc; and (v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that $\sim 60~{{\ \rm per\ cent}}$ of the proto-Milky Way galaxies are comprised by 1 dominant system (1/5 ≲M⋆/M⋆, proto-MilkyWay≲ 4/5) and 4–5 lower mass systems (M⋆/M⋆, proto-MilkyWay≲ 1/10); the other $\sim 40~{{\ \rm per\ cent}}$ are comprised by 2 dominant systems and 3–4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems.

 

The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW’s stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of the LMC's DF wake that compare the structure, kinematics, and stellar tracer response of the DM wake in cold DM (CDM), with and without self-gravity, versus fuzzy DM (FDM) withm_a= 10⁻²³eV. We conclude that the self-gravity of the DM wake cannot be ignored. Its inclusion raises the wake’s density by ∼10%, and holds the wake together over larger distances (∼50 kpc) than if self-gravity is ignored. The DM wake’s mass is comparable to the LMC’s infall mass, meaning the DM wake is a significant perturber to the dynamics of MW halo tracers. An FDM wake is more granular in structure and is ∼20% dynamically colder than a CDM wake, but with comparable density. The granularity of an FDM wake increases the stars’ kinematic response at the percent level compared to CDM, providing a possible avenue of distinguishing a CDM versus FDM wake. This underscores the need for kinematic measurements of stars in the stellar halo at distances of 70–100 kpc.

 

The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the “Timing Argument,” which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history and including the recent pericentric passage of the Large Magellanic Cloud (LMC), and recent work has found that the MW disk is moving with a lower bound “travel velocity” of ∼32 km s⁻¹with respect to the outer stellar halo. Previous Timing Argument measurements have attempted to account for this nonequilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10%–12% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume.

 

In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological-baryonic simulations. Using these potentials, we calculate actions for star particles in tidal streams around three galaxies with varying merger histories at each snapshot from 7 Gyr ago to the present day. We determine the action-space coherence preserved by each model using the Kullback–Leibler divergence to gauge the degree of clustering in actions and the relative stability of the clusters over time. We find that all models produce a clustered action space for simulations with no significant mergers. However, a massive (mass ratio prior to infall more similar than 1:8) interacting galaxy not present in the model will result in mischaracterized orbits for stars most affected by the interaction. The locations of the action space clusters (i.e., the orbits of the stream stars) are only preserved by the time-evolving model, while the time-independent models can lose significant amounts of information as soon as 0.5–1 Gyr ago, even if the system does not undergo a significant merger. Our results imply that reverse-integration of stream orbits in the MW using a fixed potential is likely to give incorrect results if integrated longer than 0.5 Gyr into the past.

 

Abstract A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N -body simulations of a spherical MW halo that recently captured a massive (1.8 × 10 11 M ⊙ ) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter particles from the present-day random distribution. Instead, the orbital poles of particles beyond R ≈ 50 kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor δ ρ max = 30% (50%) with respect to underdense regions and δ ρ iso = 15% (30%) relative to the isolated MW simulation (no LMC) between 50 and 150 kpc (150–300 kpc). The clustering appears after the LMC’s pericenter (≈50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: (1) the LMC shifts the velocity and position of the central density of the MW’s halo and disk; (2) the dark matter dynamical friction wake and collective response induced by the LMC change the kinematics of particles; (3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least ≈1:10 mass ratio), causing the clustering of orbital poles of halo tracers.